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We study by numerical simulations the two-dimensional Inelastic Maxwell Model (IMM), and 
show how the inelasticity of collisions together with the fluctuations of the number of collisions 
undergone by a particle lead to energy fluctuations that decay like a power-law. These fluctuations 
are associated to a shrinking of the available phase space. We find the asymptotic scaling of these 
energy fluctuations and show how they affect the tail of the velocity distribution during long time 
intervals. 
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I. INTRODUCTION 

Inelastic gases are low density systems composed of 
macroscopic particles themselves performing inelastic 
collisions Q, Due to their inelasticity, these sys- 
tems dissipate kinetic energy, so that their total energy 
asymptotically vanishes if they are not supplied by an 
external energy source. Nonetheless, it has been shown 
H [J H IE that their velocity distribution, when 
homogeneous, usually reach a self-similar solution, i.e. 
form preserving solution whose time dependence occurs 
through one parameter 
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T(t) is the granular temperature defined kinctically by 



(1) 



1 mV 2 



(2) 



where V=v — u is the random velocity, u is the local 
mean velocity and d is the dimension of the system. 
The average is performed with the one particle velocity 
distribution f(y;t). Such scaling velocity distributions 
have been observed in a large variety of kinetic models, 
and have been shown to generically highlight overpopu- 
lated hi gh e nergy tail, whose details depend on the model 

In Statistical Physics, interactions usually play a mix- 
ing role whose effects is to bring the system into its equi- 
librium state |ll| . This justifies, for instance in kinetic 
theory, to use the average number of collisions in order 



to measure relaxation to equilibrium. In the case of in- 
elastic gases, on the other hand, it is well known that 
the dissipative collisions prevents the system to reach its 
equilibrium state, i.e. additionally to a mixing effect that 
randomises the particles velocities, inelastic collisions act 
as an energy sink that makes their kinetic energy de- 
crease. 

In this paper, we focus on the effect of this additional 
effect on simplified kinetic models, the so-called Inelas- 
tic Maxwell Models (IMM) ©, El El El El After 
introducing the model in section IIA, we perform numer- 
ical simulations, thereby showing that the average en- 
ergy of particles depends on their collision history, i.e. 
on the number of collisions they have undergone in the 
course of time. This relation that also appears in simpler 
one-dimensional models EH i is related to the fluctuating 
character of the number of inelastic collisions. We show 
that the resulting energy fluctuations correspond to fluc- 
tuations of the phase space volume and that they relax 
like power-laws toward the asymptotic state. 



II. COLLISION NUMBER STATISTICS 

A. Inelastic Maxwell Model 

Usually, inelastic gases are defined as assemblies of 
smooth inelastic hard spheres (IHS), i.e. particles whose 
interactions do not transfer angular momentum, and are 
instantaneous and dissipative. In the absence of external 
force, the grains move freely between successive collisions 
and undergo a collision when particles i and j are in con- 
tact, with the collision rule: 
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FIG. 1: Average energy of particles as a function of their 
number of collisions, after 50 collision per particle. The 
dashed line is the constant value < E >= 1. The system 
is composed of 10000000 particles with a = 0.8, a — 0.9 and 
a = 1.0. 
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FIG. 2: Energy deviation e as a function of the rescaled de- 
viation of the number of collisions S, after 40, 70 and 100 col- 
lisions per particles (cpp). The gas is composed of 10000000 
particles with a — 0.8 and a = 0.9. 



where v.y = Vj — Vj and e is the unitary vector along 
the axis joining the centers of the two colliding spheres, 
e = t^t. The primed velocities arc the velocities before 
the collision, the unprimed ones are the post-collisional 
velocities. Energy dissipation is accounted through the 
so-called normal restitution coefficient, a £ ]0 : 1]. a—1 
corresponds to the elastic limit. Let us also stress that we 
assume that the restitution coefficient a is constant, i.e. 
we neglect the dependance of a on the relative velocity 
of the particles [Trill- 
in the low density limit, by assuming that pre- 
collisional correlations may be neglected E3 , the system 
is described by the inelastic Boltzmann equation. In the 
following, we assume that the system is and remain ho- 
mogeneous. Moreover, we use standard mean field meth- 
ods in order to simplify the mathematical structure of the 
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FIG. 3: Energy fluctuation <r e as a function of the number 
of collisions per particles. The gas is composed of 10000000 
particles with a = 0.8 and a = 0.9. 



collision operator. To do so, we assume that the collision 
frequency between particles i and j, i.e. proportional to 
Vij in the case of hard spheres, may be approximated by 
the mean field quantity y/T{t). This leads to the follow- 
ing kinetic equation: 



dt 2tt 
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dv 2 [-/( Vi ;t)/(v.;t) 
/(v <; t)/(v i5 i)] (4) 



where fi(t) = f{vi,t) and / t '(i) = /(v^t) and 6 is de- 
fined by the relation cos (9 = ^p 3 ■ Angular integrals 

and the factor \/T{t) have been absorbed into the time 
scale. This model is usually called Inelastic Maxwell 
Model (IMM), and several variations of it have been con- 
sidered in the literature [1 H H El El It is well 
known that it leads to an exactly solvable sets of equa- 
tions for the velocity moments m n {t) —< (v 2 ) n >, and 
that its asymptotic scaling solutions is characterised by 
a power-law tail , where fJ.(a) goes to infinity in 

the elastic limit a — > 1. The corresponding first diverg- 
ing moment of the scaling solution is found by solving 
a transcendental equation [l3| . Let us also remind the 

f 1 — a 2 ) 

closed temperature equation d t T — —- — that obvi- 
ously leads to an exponential decay for this quantity. 



B. Power-law relaxation and scaling 

In order to explore the energy statistics of the parti- 
cles, we have performed Direct Simulation Monte-Carlo 
(DSMC) simulations of the IMM 19]. They are based on 
the stochastic interpretation of the kinetic equation Eq0] 
and consists in picking randomly pairs of colliding parti- 
cles at each step. It is straightforward to show that the 
number of inelastic collisions N undergone by a particle 
at a given time t fluctuates. This quantity is distributed 
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according to the Poisson law , where, by construc- 
tion, t is the average number of collisions suffered by 
particles. Let us remind that the Poisson distribution 
satisfies 

< N > t = t 

/< A 2 > t -< A> 2 i 

c = W r; — 5 1 = t 2 - ( 5 ) 

Y < N > 2 w 

In order to evaluate the effect of these fluctuations on 
the particles energies, we have started simulations from a 
Maxwcll-Boltzmann initial condition, and let the simula- 
tion run during t collisions per particle. At that time, the 
simulation is paused, the particles velocities are rescaled 
so that the average energy < E >— 1 and we measure the 
average energy < E >a of particles having performed 
N = t + A collisions. The results (Fig[IJ clearly show 
that the average energy of particles is a decreasing func- 
tion of their number of collisions, thereby showing that 
ensembles of particles, discriminated by their number of 
collisions, are characterised by different quantities of en- 
ergy on average. This feature comes from the fact that 
particles having undergone less collisions have also dissi- 
pated less energy due to inelastic collisions. Simulations 
also show that this relation obviously does not take place 
when the collisions are elastic, as expected. When a < 1, 
the following scaling expression is found for long times 
t > 10 and small values of 5 = A/t (Fig|2J 



This initial increase is due to the breaking of relation 
for small times, and to the non-vanishing value of e(0). 
In the case of elastic collisions a = 1 , v< e 2 > is strictly 
zero, i.e. the dynamics are restrained on the initial con- 
stant energy surface. In contrast, when a < 1, the initial 
increase of < e 2 > corresponds to a broadening of the 
available volume in phase space, that is followed by a 
power-law shrinking of this volume. Let us stress that 
these fluctuations relax on much slower time scales than 
the usual exponential relaxations taking place in kinetic 
theory. 

Amongst others, these fluctuations contribute to the 
rapid emergence of fat tails in the velocity distribution. 
Indeed, an averaging of distributions with different mean 
energies is well-known to overpopulate the tail of the dis- 
tribution [2(], E]]. In the present case, the existence of 
these very energetic particles (as compared to the av- 
erage) arises due to the fluctuations in the number of 
collisions. We have checked this effect by DSMC. To do 
so, we have focused on particles belonging to the tail 
of the distribution, i.e. the N t particles whose energy 
is larger than k times the average, say k — 5. Among 
these Nt particles, we consider the proportion p t of par- 
ticles having performed less collisions than the average t, 
p t = Nf /Nt. The time evolution of p t shows a power- 
law, i.e. slow, relaxation that confirms the important 
role played by the observed fluctuations in the high en- 
ergy tail for very long times. 



e(A) = -r(a) j + ... (6) 

where e(A) is the deviation to the average defined by 
e =< E >a —1 and r(a) is a positive parameter that 
vanishes in the limit a — > 1. For small times, in con- 
trast, we observe a breaking of this relation, and a non- 
negligeable constant at A = 

e(A) =e(0)-e'(0)A + .... (7) 

Let us note that e(0) increases with t for small times, 
reaches a maximum around t = 5 and decreases very fast 
later on. 

Relations and [S] imply that the fluctuations of < 
E >s decrease asymptotically like power-laws 

a e = V< e 2 > ~ —r(a) i - *. (8) 

We observe this behaviour by DSMC (FigJSJ, namely 
V< e 2 > decreases like the power-law t~z after an in- 
crease that takes place during some mean collision times. 



III. CONCLUSION 

In this paper, we focus on the phase space properties 
arising in inelastic gases. To do so, we study a two- 
dimensional Inelastic Maxwell Model (IMM), which is a 
mean field approximation of the Boltzmann equation for 
inelastic hard discs. By performing DSMC simulations of 
the model, we show a statistical relation between num- 
ber of inelastic collisions undergone by a particle and 
its average energy. This relation leads to complex non- 
equilibrium features, namely the broadening of the avail- 
able volume in phase space, that is followed by a power- 
law shrinking of this volume. This mechanism is specific 
to inelastic gases, where energy is dissipated at each col- 
lision, and has no counterpart in elastic gases. Finally, 
let us stress that this work opens perspectives for a phase 
space description of inelastic gases, e.g. the use of Liou- 
ville equations for inelastic hard spheres [2^, |2^| . 
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